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Abstract: We study the extended thermodynamics, obtained by considering the cos¬ 
mological constant as a thermodynamic variable, of STU black holes in 4-dimensions in 
the fixed charge ensemble. The associated phase structure is conjectured to be dual to 
an RG-flow on the space of held theories. We find that for some charge configurations 
the phase structure resembles that of a Van der Waals gas: the system exhibits a family 
of first order phase transitions ending in a second order phase transition at a critical 
temperature. We calculate the holographic entanglement entropy for several charge con¬ 
figurations and show that for the cases where the gravity background exhibits Van der 
Waals behavior, the entanglement entropy presents a transition at the same critical tem¬ 
perature. To further characterize the phase transition we calculate appropriate critical 
exponents and show that they coincide. Thus, the entanglement entropy successfully 
captures the information of the extended phase structure. Finally, we discuss the phys¬ 
ical interpretation of the extended space in terms of the boundary QFT and construct 
various holographic heat engines dual to STU black holes. 
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1 Introduction 

The gauge-gravity correspondence [1-3] has uncovered deep connections between black 
holes in anti-de Sitter space and quantum field theories. According to the AdS/CFT 
dictionary, these solutions are the holographic duals to strongly-coupled large-A gauge 
theories at finite temperature. Since black holes have thermodynamical properties, they 
seem to undergo phase transitions in the same manner as ordinary thermodynamical 
systems, and the implications of such phase transition for the boundary theory are worth 
exploring. The best known example of such a transition is the Hawking-Page transition 
[4], a first order transition between the Schwarzschild-AdS solution and thermal AdS. 
The Hawking-Page transition was interpreted in the holographic context by Witten [5] 
as the gravity dual to the confinement/deconfinement transition of AT = 4 SYM. 

The study of black hole transitions was then extended to include charged black 
holes in AdS [6, 7]. Compared to the Schwarzschild-AdS case, the phase structure of the 
Reissner-Nordstrom-AdS solution was found to be richer and dependent on the chosen 
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statistical ensemble. In particular, in the canonical ensemble (obtained by keeping fixed 
the electric charge of the black hole), the behavior of the inverse temperature versus 
the horizon radius (or, equivalently, the black hole entropy) reveals a phase transition 
similar to that of a Van der Waals gas. Moreover, there is a critical value of the charge 
where the system has a second-order critical point. 

In a recent series of papers, the analogy with the Van der Waals gas has been made 
more precise in the framework of extended black hole thermodynamics [8-25]. In this 
context, the cosmological constant is treated as a dynamical variable and is associated 
with the pressure of the system, P = —A/8ttGn, while its conjugate quantity is identified 
as the thermodynamical volume V. The first law is then extended to 1 


dH = TdS + VdP , 


( 1 . 1 ) 


where P[ is the enthalpy of the system (or, equivalently, the ADM mass M), T = /3 1 is 
the black hole temperature and S is the Bekenstein-Hawking entropy, 


S 


A 
4 Gjy 


( 1 . 2 ) 


The extended phase space and associated thermodynamics have been studied for the 
RN-AdS solution, Taub-NUT-AdS, Taub-Bolt-AdS [26-28] and their cousins in higher 
derivative theories such as Gauss-Bonnet and Lovelock gravity [29-34], In simple cases 
such as the RN-AdS solution, the thermodynamical volume is shown to coincide with a 
naive integration over the black hole interior (in the Schwarzschild slicing), 


Gp- 2 r+ 1 

D - 1 


(1.3) 


However, it remains unclear how should it be interpreted, since the volume of a black 
hole is not a well-defined quantity and depends on the slicing of the spacetime. 2 Nev¬ 
ertheless, in most of the cases studied it was found that isotherms on the black hole 
PV diagrams exhibit the characteristic Van der Waals behavior: the system undergoes 
a family of first order phase transitions ending in a second order phase transition at a 
critical temperature. 

Unlike the Hawking-Page transition, the Van der Waals transition unfortunately 
lacks a clear field theory interpretation. The first steps in this direction were made 
recently in [35]. In this paper, the author proposed that varying the cosmological con¬ 
stant in the bulk corresponds to perturbing the dual CFT, triggering a renormalization 
group flow. The transition is then interpreted not as a thermodynamical transition but, 

x For charged black holes there is an extra term that must also be included into the first law, dH = 
TdS + VdP + QdQ. However, in the canonical ensemble dQ = 0 so it reduces to (1.1). 

2 Intriguingly, the thermodynamic volume can be negative in the case of the AdS-Taub-NUT space- 
time. In [26, 27] the volume variable was interpreted in terms of the formation process of the spacetime. 
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instead, as transition in the space of field theories. This statement can be intuitively 
understood since, in holographic theories, the cosmological constant is usually associ¬ 
ated to the rank of the gauge group of the dual theory N , or its central charge c. Thus, 
changing the value of A is equivalent to changing the boundary theory. Nonetheless, it is 
worth emphasizing that the details of the field theory interpretation are still, to a large 
extent, open questions. In particular, the precise role of the thermodynamical volume 
V remains to be fully understood. 

The Van der Waals transition has also been observed in a different context [36], 
although no connection with extended thermodynamics was pointed out. The authors 
showed that, like black hole entropy, entanglement entropy also undergoes a Van der 
Waals-like transition as the temperature varies. However, the cosmological constant 
was kept fixed, so the emphasis was in finding potential holographic applications to 
condensed matter systems, e.g. [37-39]. Here we point out that there is indeed a 
connection between these two approaches. In particular, we show that the relevant 
parameter that determines the transition is a dimensionless combination between the 
temperature and the pressure, ( = /3 2 P, which can be tuned either by varying the 
cosmological constant (changing the theory) or the temperature (changing the state). 
Since the transition is also visible in the entanglement entropy computation, this result 
suggests that we can use this field theory observable as an efficient tool to unravel the 
full space of theories spanned in the PV space. 

Before proceeding further, let us motivate the choice of gravity solutions that we will 
consider throughout this paper. On the gravity side, the idea of varying the cosmological 
constant finds its natural setting in theories where the cosmological constant arises as 
a vacuum expectation value. For example in stringy black hole solutions where the 
cosmological constant is the minimum of the potential of scalar fields. This was the 
main motivation for [40]. In that work, the thermodynamical volume of a class of 
supergravity solutions called STU black holes was computed and studied in detail. The 
authors found that for STU black holes the thermodynamical volume does not coincide 
with a naive integration over the black hole interior. 3 Another advantage of studying 
the extended phase structure in a specific top-down constructions such as the STU black 
holes is that the held theory interpretation is easier to understand. In particular, the 
scalar fields that appear in the bulk system are dual to certain operators in the dual 
theory that are turned on and off and, therefore, triggering a nontrivial RG-how in the 
boundary. Thus, in this work, we will focus on the STU black hole solutions. Other 
properties of the STU black holes have been explored in [41-44], 

In the present work we study the extended thermodynamics of STU black holes. 
We first ask ouserlves how generic is the Wan der Waals behavior in these more general 
backgrounds that involve four charges. We find that it is not quite generic; it is not 

3 Rather, the thermodynamical volume agrees with the integral of the scalar potential over the black 
hole interior [40]. 
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present for all charge configurations. It is then natural to ask whether the holographic 
entanglement entropy correctly captures this property. The answer is yes. We show 
that -even for small entangling regions- the holographic entanglement entropy not only 
distinguishes between STU black holes with and without Van der Waals but also exhibits 
the same critical temperture and the same critical exponents as the ones obtained in the 
extended thermodynamic framework. 

The paper is organized as follows: in Section 2, we start by briefly reviewing the 
thermodynamics and phase transitions of STU black holes for different charge configu¬ 
rations. We focus on the case of 4-dimensional black holes, although we suspect that the 
conclusions must also apply to similar spacetimes in 5 and 7 dimensions. In Section 2.2, 
we study the isotherms on the PV diagram and show that the Van der Waals behavior 
is observed in two special configurations: the cases of three and four equal charges. In 
Section 2.3 we study further the phase transitions and compute the characteristic critical 
exponents of various thermodynamical variables around the critical point. In Section 
2.4, we comment and clarify some previous ideas on the field theory interpretation of 
the PV phase space and the construction of holographic heat engines. In Section 3, we 
turn to the computation of holographic entanglement entropy. We verify that for charge 
configurations that exhibit a Van der Waals behavior, the entanglement entropy indeed 
undergoes a similar transition at the same critical temperature. In Section 3.3 we com¬ 
pute the critical exponent of the entanglement entropy around the phase transition and 
show that it is in agreement with the expected result for a Van der Waals gas. Finally, 
we close in Section 4 with conclusions and point some possible directions for future work. 

2 Extended thermodynamics of STU black holes 

2.1 Background and thermodynamical variables 

The STU black hole in 3+1 dimension is a black hole solution of gauged N = 8, D = 4 
supergravity, first derived in [42, 43]. Gauged N = 8, D = 4 supergravity arises by the 
compactihcation of D — 11 supergravity on S 7 , then truncating to the massless sector. 
The gauge group, which is also the isometry group of the internal space, is SO( 8). In 
order to find black hole solutions, the gauge group is further truncated to its Cartan 
subgroup 17(1) x 17(1) x 17(1) x 17(1). Hence, the black hole solution has four different 
electric charges. 4 The bosonic part of the 4-dimensional effective action reads: 

s = ^j ^7+ - \ E U+T - v - \ y x 7 xt) . (2-1) 

4 Analogous solutions exist in 4+1 and 6+1 dimensions. The 4+1 solution was obtained from com- 
pactifying D = 10 Type IIB supergravity on S 5 . Since the isometry group of S 5 is SO(6), with Cartan 
subgroup 17(1) x 17(1) x 17(1), the 4+1 black hole has three electric charges. The 6+1 solution was 
obtained from compactifying D = 11 supergravity on S 4 . Since the isometry group is £>0(5) and its 

Cartan subgroup is 17(1) x 17(1), the 6+1 black hole has two electric charges 
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where 


X! = e~^, 


( 2 . 2 ) 

CZi ( 1 , 1 , 1 ) , 


(2.3) 

Pi 

to 

1 

1— 1 

1 


(2.4) 

®3 ( 1 ) I; 1 ) j 


(2.5) 

O 4 ( 1 , 1 , 1 ) , 


( 2 . 6 ) 

and the scalar potential is given by 



v= -g 2 J2 x ‘ x i- 


(2.7) 

i<j 



The solution discovered in [43] is given by: 



ds 2 = -(H^HsH^-^fdt 2 + (tfitf 2 tf 3 tf 4) 1/2 | 


( 2 . 8 ) 

A < = V q i(* + *rn) dt 
r + qi 


(2.9) 

v 

X '~ H, 

where 

2 m 0 0 

/ = 1- + g 2 r 2 H 1 H 2 H 3 H A , 

r 


( 2 . 10 ) 


( 2 . 11 ) 

Hi = 1 + Q - . 

r 


( 2 . 12 ) 

The relevant thermodynamic quantities are given by: 



H = m+ \ ^2 qi > 
i 


(2.13) 

Qi = + 2m ) > 


(2.14) 

^ ~ n ^ r ++* > 

i 


(2.15) 

r = ^n^ 1 / 2 (r + ), 

i 


(2.16) 

_ yJqMi + Zm) 

2 (r+ + ft) 


(2.17) 

3 2 

P = —g , 

8 tt ’ 


(2.18) 


5 
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(2.19) 



where H is the enthalpy (identified with the ADM mass of the spacetime), Qi are the 


physical charges, S is the entropy, T is the Hawking temperature, <3>j is the electrostatic 
potential at infinity, P is the pressure and V is the thermodynamic volume. Notice that 
we are working in Planck units (c = H = ks = G = 1), so all quantities are dimensionless. 
The outer horizon is located at r = r + , the largest root of f{r + ) = 0. 

For the STU black hole, the first law of thermodynamics reads: 
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( 2 . 20 ) 


2=1 


In the remainder of the paper we will work in the canonical ensemble, i.e. we keep all 
the physical charges Qi fixed, so the first law simplifies to: 


( 2 . 21 ) 


dH = TdS + VdP. 


We can easily check that the Hawking temperature (2.16) and the thermodynamic vol¬ 


ume (2.19) can be recovered from the enthalpy by the usual thermodynamic relations: 



( 2 . 22 ) 



(2.23) 


2.2 Extended phase structure 

We will now proceed to study some aspects of the extended phase structure of STU 
black holes. In order to fix notation and to illustrate the procedure we will first review 
the special case of AdS-RN which have already been studied in various papers, e.g. [35]. 

2.2.1 The AdS-RN black hole 

An AdS-RN black hole is a special case of an STU black hole when all the charges are 
equal. Setting q 1 = q 2 = q?, = = q in (2.8) we obtain, 



(2.24) 


where 



(2.25) 
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The thermodynamic quantities entering (2.21) are now, 


H — m + q , 



S = n(q + r + ) 2 , 

4 

V = -7r (g + r+) 3 , 


mr + — g 2 (q — r + )(g + r + ) 3 


27rr + (g + r + ) 2 


(2.26) 

(2.27) 

(2.28) 
(2.29) 


(2.30) 


A peculiarity of the AdS-RN solution is that the pressure and volume determine 
each other. Indeed, from equations (2.28) and (2.29) it follows that: 



(2.31) 


One might be worried about this, since it means that the volume and entropy cannot 


be used as independent thermodynamical variables. Fortunately this entropy-volume 


degeneracy is specific to AdS-RN, volume and entropy are independent variables in the 
other STU special cases (as well as for rotating black holes, for instance). Also recall 
that the physical charge is not q but 



(2.32) 


At this point would be convenient to introduce the radial Schwarzschild coordinate 


R = r + q which would allow us to write the metric in a more familiar form 5 and 


obtain the equation of state in a straight forward manner. However, in the general case 
such strategy is not available and one has to proceed numerically in order to obtain the 
equation of state. Therefore, we will refrain here from switching to the Schwarzschild 
coordinate here and work directly with equations (2.24)-(2.30). From (2.29) and (2.25) 
we obtain horizon and mass parameter as a function of q, V and P: 



(2.33) 


(2.34) 


5 In these coordinates the AdS-RN metric takes the form 


ds 2 = -F(R)dt 2 + + R 2 dP 2 , 

t [K) 


with 
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0.2 0.4 


Figure 1. PVT diagram for AdS-RN with charge Q = 0.05. The solid color curves are 
isotherms. The blue one clearly exhibits Van der Waals behavior but the green one does 
not. The critical temperature at which the transition occurs is depicted in red, and can be 
determined from the conditions (2.36). 


Using (2.34), (2.30) and (2.32) we can then solve for the charge parameter q as a function 
of the physical charge Q , the volume V and the pressure P. 

With all these ingredients at hand, the temperature (2.27) defines the equation of 
state, T = J r (Q , V, P), which we can rewrite as 

P = A UWV—_ m 2/3 J_ + *91) (2 35) 

87T y 3 / \4n) V 2 / 3 + V 4 /3 J ' 


As shown in Figure 1, this equation of state exhibits a Van der Waals-like phase transition 
at sufficiently low temperature. From this plot, we can see that there is a critical 
temperature, T c , at which there is an inflection point with 

/ dP\ _(d 2 P\ 

\dv) T -\dv 2 ) T - 

This corresponds to a second order phase transition, above which there is no longer a 
Van der Waals transition. The critical pressure, volume and temperature are given by 

Pc = oo) 1 ^.o ’ Vc = 64V67rQ 3 and /3 C = qVQhQ , (2.37) 

3o47rqV 

respectively. 

Before proceeding further let us note that it may appear surprising that the dual held 
theory, which lives on a finite volume, could undergo a phase transition at all. Phase 
transitions require non-analyticities in the partition function. Such non-analyticities 
only arise in system with infinitely many degrees of freedom, which, in more familiar 
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(2.36) 







thermodynamical systems, implies infinite volume. In our case, there are indeed an infi¬ 
nite number of particles since we are working in the large -N limit, hence the possibility 
of phase transitions. 

Notice that Van der Waals behavior of AdS-RN was first discovered [ 6 ] without 
considering the cosmological constant as a thermodynamical variable. Indeed, if we do 
not consider the extended thermodynamics and work at a fixed pressure we arrive to 
the following equation of state: 


T = 


1 

47T 



(2.38) 


It is possible to show that there exists a critical value of the charge Q c at which this 
curve has an inflection point 


= (&l\ 

dSj Qc \dS*) Qi 


(2.39) 


for which one finds 


Sc 


7r 

6 ?’ 



and 



(2.40) 


For Q > Q c , the curve is thermodynamically stable, while for Q < Q c , the plot develops 
an unstable branch in addition to two stable branches. Also, note that if we plug g(Q c ) 
in (3 C we recover the same critical temperature as in (2.37). This indicates that the 
relevant (dimensionless) parameter that determines the transition is £ = (3g = (TL) -1 , 
which can be tuned either by varying the AdS radius L (i.e. the cosmological constant) 
or the temperature T. 


2.2.2 Other charge configurations 

Our goal here is to repeat the analysis outlined in the previous subsection for different 
charge configurations. The logic and the steps are the same but once we move away 
from the AdS-RN case it is not possible to find an analytic expression for the equation 
of state. Thus, we have to proceed numerically. 

We will focus on cases where one, two or three charges are turned on but, for the 
sake of simplicity, we will consider all of them to be equal, i.e. we will set qi = q Vi. The 
thermodynamic quantitites are given in (2.14)-(2.19). In particular, the temperature for 
a general STU black hole is (2.16), 


T = 


2 P 


3 A /n: =1 (i+j) 


n 

i¥=i 


i —1 


i 

r + 


+ 


2 m 3 


87 tP 


+ 2r +n 


2—1 


i+® 

r+ 


(2.41) 
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Figure 2. Isotherms in the PV diagrams of the STU black hole in 3+1 dimensions. Top left: 
The single charge case Q\ = 0.05 and Qo = Q 3 = Qa = 0. The temperatures chosen are T = 
0.55 (green), 0.5 (orange) and 0.45 (blue). Top right: The two-charge case Q\ = Q 2 = 0.05, 
Q 3 = Qi = 0- The temperatures chosen are T = 0.96 (green), 0.93 (orange) and 0.90 (blue). 
Bottom left: The three-charge case Qi = Q 2 = Q 3 = 0.05, Q 4 = 0. The temperatures chosen 
are T = 0.56 (green), 0.5391 (orange) and T = 0.51 (blue). Bottom right: The AdS-RN case 
with Q 1 = Q 2 = Q 3 = Q 4 = 0.05. The temperatures chosen are T = 0.48 (green), 0.43316 
(orange) and 0.4 (blue). For the two-charge, three-charge and AdS-RN cases, the orange curve 
is at critical temperature. 


Using (2.14), (2.19) and (2.11) one can solve for r + , m and qi in terms of Qi and V. 
Thus, (2.41) implicitly defines the equation of state relating T, P, V and the physical 
charge Q. In all cases studied it is possible to obtain three dimensional PVT plots similar 
to Figure 1. However, for maximum clarity we choose to plot constant temperature 
curves in a PV diagram. 

The numerical results are shown in Figure 2 , for different number of charges. As 
mentioned before, we only consider the cases where the nonzero charges are all equal 
to each other. However, our numerical explorations indicate that the behavior is qual¬ 
itatively similar if the nonzero charges are not equal to each other. That is to say, the 
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AdS-RN case will be representative of all cases where all four charges are nonzero, etc. 
As can be seen from the plots, the Van der Waals behavior is present in the case of 
three and four active charges, i.e. below a certain critical temperature T c the isotherm 
has two stable branches (where the volume decreases with increasing pressure) and one 
unstable branch (where the volume increases with increasing pressure). Moreover, for 
the three-charge case we can use the 3-dimensional plot P(V,T ) to find T c , P c and V c 
numerically. For Q = 0.05, the results are: 


(T c , V c , P c ) = (0.539, 0.0264, 0.524). 


(2.42) 


But from dimensional analysis, we know that: 


1 o 1 

T c oc — , V c oc Q J , P c oc — , 


(2.43) 


where the proportionality factors are dimensionless numbers. Therefore, we can easily 
deduce the critical temperature, volume and pressure for arbitrary Q: 


T r 


0.0270 

Q 


V c = 211.254 Q 3 


and 


Pr 


0.00131 

~Q^~ 


(2.44) 


respectively. These are equivalent expressions to (2.37), for the case of the AdS-RN 
black hole. Rs worth emphasizing that if we change the value of some of the charges, 
the constants of proportionality will change, but the phase structure will be similar. 

The Van der Waals behavior is absent in the two other cases (one or two active 
charges). In the case of two active charges, at sufficiently low temperatures, the isotherms 
has an unstable branch at small volumes and a stable branch at large volumes. As 
the temperature increases, the system undergoes a phase transition at some critical 
temperature where the unstable branch is squeezed out. The critical temperature at 
which this happens is T c = 0.0465Q -1 . Finally, in the case of only one active charge, 
there is an unstable branch for small volume and a stable branch at large volumes, and 
the unstable branch is present at all temperatures. This is the same behavior as for 
the zero-charge case (i.e. Schwarzschild-AdS), so there must also be a transition at 
small volumes reminiscent of a Hawking-Page transition, presumably to a fixed charge 
ensemble in pure AdS. 

Next, we turn our attention to plots of the temperature versus the black hole entropy. 
In Figure 3, we present these plots for the four different cases of the STU black hole that 
we are considering. Each curve corresponds to a different value of the charge, while keep¬ 
ing the AdS radius fixed L = 10. As previously, we have generated 3-dimensional plots 
of (S, (3, Q) but we have chosen to present 2-dimensional slices for clarity. Remarkably, 
the plots are qualitatively the same as Figure 2: for the three-charge and four-charge 
cases, the curves have an unstable middle branch at sufficiently low temperature. 6 We 

6 This branch is indeed unstable, since its specific heat Cp = T (§§) p is negative. 
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Figure 3. Plots of the inverse temperature as function of the black hole entropy S for the one- 
charge case (top left), two-charge case (top right), three-charge case (bottom left) and AdS-RN 
case (bottom right). For each panel, each of the 3 curves correspond to a fixed numerical value 
of the common charge. In the upper two plots, the middle curve corresponds to the critical 
charge and the red horizontal line indicates the critical temperature. 


indicate the critical temperature in each case by a red horizontal line. For the AdS- 
RN case, the critical charge and critical temperature follow from (2.40): Q c = 5/6 and 
Pc = 10 a / 3 / 27 t . For the three-equal-charge case, we made use of the 3-dimensional plot 
to solve for the critical temperature, entropy and charge: 


(Q c , Pc, S c ) = (1.048, 38.871,46.959). (2.45) 

These values are valid for g = 1/L = 0.1. From dimensional analysis, we know that: 


Qc oc - , 
9 


Pc<x~ , 
9 


S c oc — . 

r 


(2.46) 


Therefore, we can deduce the formula for the critical charge, temperature and entropy 
at arbitrary pressure: 


Qc = 


0.1048 

9 


Pc = 


3.8871 

9 


and 


S c = 


0.46959 
9 2 


(2.47) 
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respectively. On the other hand, the two-charge and single-charge cases do not exhibit 
the Van der Waals transition. The two-charge case also undergoes a phase transition 
where the unstable branch is shrunk to zero. This happens at the critical charge Q c ~ 
0.20704 /g and critical temperature f3 c ~ 4.44187/g. Finally, for the one-charge case, the 
unstable branch persists at all temperatures. 

It is worth pointing out that for charge configurations that exhibit a transition, 
the parameter £ = /3g = (TL) _1 is fixed in both scenarios, i.e. in the PV diagram at 
constant T or in the TS diagram at constant P. This result shows that we can study 
the phase transition either by tuning the pressure or the temperature, with completely 
different interpretations in terms of the dual field theory. We will comment more on this 
point in Section 2.4. 


2.3 Characterization of the phase transitions 

An important characterization of continuous phase transitions is given by the so-called 
critical exponents, which determine the leading power dependence of various quantities 
close to the phase transition. Remarkably, thermodynamical systems are often found 
with identical critical exponents even though their microscopic details may differ. This 
suggests that the critical exponents may be used to classify phase transitions into univer¬ 
sality classes. In particular, from the equation describing the critical isotherm P(T C , V ), 
we can extract the critical exponent 5 defined to be: 

\P-Pc\ oc \V- VI 5 . (2.48) 


The factor of proportionality in the equation above depends on the microscopic details 
of the system, and therefore is unimportant. For the AdS-RN case equation (2.48) takes 
the form 

\P(T C , V) — P c \ = -- |V-V C | 3 , (2.49) 

' ' 2916v / 67t 4 (4Q) 11 ' 1 V ' 


and 5 is found to be [11]: 


5 = 3 


(2.50) 


Next, we consider the specific heat at constant pressure Cp, defined by: 


C P = T 



(2.51) 


We will denote a as the critical exponent associated with Cp, which at leading order 
behaves as 

C P oc (T - T c )~ a . (2.52) 

For AdS-RN, we have the expansion: 


1 / 3 ( 5 , Qc) ~ Pel = 192^/t rP 5 / 2 (S - S c ) 3 , 


(2.53) 



and after some algebra, we find that the numerical value of a is given by [36] 

a= \- ( 2 - 54 ) 

Note that, since we are working in the canonical ensemble, the specific heat we refer to 
as Cp is also at constant charge Q. More generally, for generically chosen paths Q(T) 
approaching the critical point, the critical exponent associated to the specific heat Cp 
its found to be 7 = 1 [11]. It is only for the special case we consider here, Q = Q c , 
that the exponent is given by 2/3. Remarkably enough, these two critical exponents 
happen to agree with the expected values for a Van der Waals gas . 7 In addition, several 
other critical exponents for AdS-RN as well as for most AdS black holes also agree with 
expectations from mean-field theory [11], The idea here is to perform a similar study 
for the STU black holes. 

Recall that for other charge configurations we do not have analytical expressions for 
the equation of state so we have to proceed numerically. We will focus on the three- 
equal-charge black hole, which is the only other case in consideration that exhibits a 
Van der Waals-like transition —see Figure 2. In order to extract the value of S we first 
take the logarithm of both sides of (2.48), 

log|P-P c | =5\og\V- V c \+C, (2.55) 

where C is an additive constant. Therefore, we can plot log |P — P c | as a function of 
log|V — V c \ in the vicinity of V c and as a result we expect a straight line with slope 
5. We present this plot in Figure 4(a), with the common value of the charges set to 
Q = 0.05. A linear fit of the data points yields: 

log |P - P c | = 2.9476 log | V - V c | + 7.0200 , (2.56) 

with a slope that suggests the exact value: 

5 = 3. (2.57) 

Notice that this is the same value that we found for AdS-RN. For comparison, the analog 
of equation (2.56) for AdS-RN with four all charges set to Q — 0.05 reads: 

log |P - P c | = 3 log | V - V e \ + 4.25105 . (2.58) 

The ^-intercept of the equation above differs from that of (2.56), even though they are 
of the same order of magnitude. This difference reflects the different microscopic details 
between AdS-RN and the three-equal-charge STU black hole. 

'Our a is not to be confused with the critical exponent usually denoted by a in statistical mechanics. 
This latter is associated with Cy. 
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Figure 4. Left panel: log |P — P c \ vs. log \V — V c \ for the three-equal-charge black hole with 
Q = 0.05. The linear fit is based on 10 points closest to the inflection point, in the volume 
range V G [0.0210,0.0254]. The critical volume is at V c = 0.0264. Right panel: log |/3 — (3 C \ vs. 
log \S — S c \ for the three-equal-charge black hole with Q = 0.05. The linear fit is based on 23 
points closest to the inflection point, in the entropy range S G [50.4622,55.3306]. 


Next, we can repeat the same procedure for the TS diagrams, and consider the 
power dependence of S around the critical point (on the critical charge curve). To do 
so, we plot log |/3 — j3 c \ vs. log \S — S'd in Figure 4(b), with g — 0.1. A linear fit of this 
curve yields: 

log |/3 - /3 C | = 3.0486 log | S- S c \ - 11.2609 , (2.59) 

suggesting that 

(3 —/3 C (x \S — S c \ 3 . (2.60) 

Then, from (2.52) it follows that 

<*=1- (2.61) 

Again we find the same value as the AdS-RN case. For comparison, the analog expression 
to (2.59) for AdS-RN, with g = 0.1, is: 

log |/3 — /3 C | = 3 log [S' — S'cl - 10.997. (2.62) 

2.4 Field theory interpretation 

The analysis of [35] lead to some puzzles regarding the implications of the extended 
thermodynamics in the context of the AdS/CFT correspondence. Here we will review 
the main issues and discuss the resolutions that have been proposed in terms of a specific 
top-to-bottom construction, where the physical interpretation is neater. 

According to the holographic dictionary, (3 + l)-dimensional asymptotically AdS 
black holes are dual to CFTs at finite temperature in 2 + 1 dimensions. However, the 
pressure and volume appearing in the extended thermodynamics are not identified with 
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the pressure and volume of the dual field theory, instead, they are interpreted as (non¬ 
thermodynamic) dynamical variables in the space of field theories [35]. 

Let us first recall the main entries in the path integral formulation of AdS/CFT. 
For a fixed cosmological constant A, the Euclidean path integral in the gravity side, 

I E /P = -logZ/P, (2.63) 

is equivalent to the Helmholtz free energy F = U — TS of the dual CFT. However, in 
top-to-bottom constructions (as the one we consider here), the value of A is set by the 
potential V(</>*) where the </>* are dynamical fields in the gravity side. These fields can be 
excited by turning on and off the charges of the system, thus, allowing the cosmological 
constant to vary. In this case, the natural quantity that we should identify with (2.63) 
should be the Gibbs free energy, G = U — TS — PV since, from the point of view of the 
extended thermodynamics, varying the cosmological constant is equivalent to varying 
the pressure P = —A/8ttGn] here A = —3/L 2 and Gn is the Newton’s constant. 

The STU black holes considered here are solutions of gauged N = 8, D = A super¬ 
gravity which arises as a truncation of D = 11 supergravity compactihed on a S 7 . This 
theory can be uplifted to full M theory and hence, these black holes can be thought of 
as arising from different configurations of M2-branes. In the effective description, the 
value of L is set by the value of the Planck length in the eleven dimensional theory Ip, 
and the number of branes N. The worldvolume theory is described in terms of a gauge 
theory with symmetries specified by the specific brane configuration; typically N is the 
rank of the gauge group so it determines the number of degrees of freedom the theory. 
Moreover, the Newton’s constant Gn also depends on N. In fonr dimensions, these two 
quantities scale as L ~ N 1 ^ and Gn ~ N~ 7 ^ [1] so one finds that the pressure depends 
nontrivially in N, specifically as P ~ N 5 ^ 6 . Is in this sense that, varying the cosmolog¬ 
ical constant A (and hence the pressure), is equivalent to changing the field theory to 
which the bulk background is dual to. 

One might wonder how this notion of RG-flow fits in the standard picture of the 
AdS/CFT correspondence. According to the UV/IR connection, the radial direction in 
the bulk maps to an energy scale in the boundary [45-48] . This emergent direction can be 
thought of as way to encode the RG-flow in a geometrical way, since it links the UV and 
the IR of the theory (which correspond to the near-boundary and near-horizon portions 
of the geometry, respectively). However, in the framework of extended thermodynamics 
is the cosmological constant itself that plays the role of an energy scale, given that its 
value fixes the AdS radius L. We argue that the physical interpretation in each case 
is different. In order to understand better the distinction of these two scenarios let us 
review other relevant entries of the holographic dictionary. One of the main properties 
of the AdS/CFT correspondence is the map between bulk fields <f)(x,r ) and operators 
in the gauge theory side 0(x). According to [2, 3], we can identify 

ZU o] = ( e -fdGMPO(x)\ (2.64) 

\ / CFT 
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where the term in the left hand side is the string partition function with boundary 
condition 8 </> —» r A ^ 3 0 o as r —» oo and the term in the right is the generating functional 
of correlation functions of the gauge theory. Notice that, for the right-hand side of (2.64) 
to be well-defined, O has to have conformal weight A, so that 0(Xx) = A ~ a O(x). More 
in general, near the boundary, solutions to the equations of motion take the general form 

A B . 

+ ( 2 - 65 ) 

where A and B are identified with the source and expectation value of the dual operator 
O, respectively. For the STU background, expanding the scalar fields (2.2)-(2.10) around 
r —> oo lead to 


0i 

02 

03 


gi + ?2 - -ql -qi + gj + ql 

2 r 4r 2 

gi ~ ?2 + g 3 - & - 9 i + ql - g 3 + ql , 

2 r 4r 2 

gi ~ ?2 - g 3 + ?4 -ql + ql + qj-ql 

2 r 4r 2 


( 2 . 66 ) 


(2.67) 


The conformal dimension of the three dual operators C?j, for i = {1, 2, 3}, is found to be 
A = 2. 9 To see this, we can alternatively expand the potential (2.7) as 

V = -6 g 2 - s 2 (0? + 0^ + 01) + C?(0?), (2.68) 


from which we can read off the AdS radius L — 1/g and the three masses m 2 = —2g 2 . 
Then, we use the formula 


A 


d + Vd 2 + Am 2 L 2 

2 


2 . 


(2.69) 


Thus, by giving different values to the charges in the bulk geometry we are, at the same 
time, turning on and off the couplings A,: of three relevant operators in the gauge theory. 
Schematically, the Lagrangian of the dual theory is modified according to 


3 

i= 1 


(2.70) 


where £o ~ .9 ymT r [F 2 + • • • ] is the Lagrangian of the original theory (with all the 
charges/scalars turned off). Furthermore, regardless of the charges, the bulk geometry 
is asymptotically AdS, and thus the physics in the UV is still governed by the CFT 

s Here assume the standard quantization. 

9 In the alternative quantization the normalizable and non-normalizable modes interchange roles so 
A = 1. Either way A < 3 and the dual operator is relevant. 
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prescribed by Cq. However, the relevant operators O, can trigger an RG-flow to a new IR 
fixed point governed by a different theory. This effect is encoded in the radial dependence 
of the bulk geometry. Indeed, one of the effects of the non-trivial profile of the scalars is 
that the bulk solution is no longer invariant under the scale transformations x M — > Xx 
r —» A -1 r and therefore the solutions might have a completely different behavior deep 
in the bulk interior. 

Now, we have argued that varying the cosmological constant provides an alternative 
notion of changing an energy scale of the dual theory. Can this notion be equivalent 
or, at the very least, related to the standard RG-flow encoded by the radial direction 
in the bulk? We believe that the answer is no. Rather, it provides a different and 
orthogonal direction in which we can move from theory to theory. Consider for example 
the AdS-RN black hole. In this case all the scalars are turned off so the Lagrangian of 
the theory is just Co, with no additional deformations. Here, it is clear that the non¬ 
trivial dependence on the radial coordinate just reflects the fact that we are not in the 
vacuum of the theory, but instead in an excited state with additional energy scales such 
as temperature and/or charge density. However, we can still think of varying the rank 
of the gauge group N, which can be achieved by changing the value of the cosmological 
constant A (or, the pressure P ). Hence, we can still have a non-trivial RG-flow. One 
issue that remains unclear in this context is how the volume V should be interpreted in 
terms of the dual theory. In the case in which we have a non-trivial RG-flow along the 
radial direction, the relevant deformations O t break explicitly the conformal invariance 
and it is natural to define the theory at a cutoff scale A c . 1(l The dependence of physical 
quantities such as the couplings A*(A C ) on the cutoff is encoded on the radial profile 
of the bulk solutions. Similarly, if we vary N (or the pressure P) there should be an 
analogous scale A* that determines how the physical N(A*) runs along the RG-flow. 
This is precisely encoded in the P(V) diagram at fixed T. Therefore, we can think of 
the volume V as a cutoff scale that determines running of N with energy. 

Schematically, then, the PV space can be thought of as a space theories with different 
N and A*; one parameter takes us from theory to theory while the other moves us along 
the RG-flow. We can, for instance, think of constructing “thermodynamic” cycles in 
this space of theories. These processes can be thought of as holographic heat engines, 
in the same sense of [35]. In Appendix A we show details of such a construction for the 
particular case of STU black holes. 

Before closing this section let us comment on the nature of the phase transitions 
in the framework of the extended thermodynamics. We have seen that for the several 
cases studied, the relevant parameter that determines the transitions is the dimensionless 
parameter ( = /3 2 P. This provides us with an interesting possibility, namely, we can 
study the transitions by either by tuning the pressure (changing the theory) or the 

10 If all the Xi vanish as in the case of AdS-RN, conformal invariance is still broken due to the presence 
of temperature and/or charge density. In this scenario it is still sensible to define a cutoff scale A c . 
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temperature (changing the state). In terms of the TS diagram at fixed P, the transitions 
occur to guarantee thermodynamic stability. In particular, the unstable branches with 
positive slope must be replaced by straight lines in order to satisfy the positivity of 
specific heat at constant pressure, 



(2.71) 


Notice that this interpretation is valid for both, the bulk and the boundary theories. On 
the other hand, the analogous statement for the case of the PV diagram at fixed T is 
given in terms of the positivity of the isothermal compressibility 11 



(2.72) 


where in this case the inequality stands for mechanical stability. The unstable branches 
in the PV diagram should replaced by a straight line obeying the Maxwell’s equal-area- 


law. However, since P and V are not related to the physical pressure and volume in the 
boundary theory, this interpretation does not hold in the dual QFT. Here, the analogous 
criterium for stability should be related to the monotonicity of N(A*) along the RG- 
flow! Exploring the consequences of this statement is clearly an interesting topic but lies 


beyond the scope of this paper. We will study it in more detail in the near future [49]. 


3 Holographic entanglement entropy 

Having established that, depending on the charges, there are different kind of transi¬ 
tions in the extended phase space of STU black holes, we want to investigate how such 
transitions are reflected in the dual held theory. As mentioned in the introduction, the 
work of [36] lead to the conclusion that at least for the AdS-RN case a similar transition 
to the one observed in the TS diagram is present in the holographic computation of 
entanglement entropy. However, in that paper no connection with the extended thermo¬ 
dynamics was made. In this section, we will show that this transition is indeed the same 
as the transition appearing in the PV space, implying that this held theory observable 
can be used as a tool to diagnose the extended phase structure of the theory. We will 
perform the calculation for various STU black holes and show that the same result holds 
more generally, thus making the conclusion more robust. 

3.1 Review of holographic entanglement entropy 

We will begin with a brief review of the Ryu-Takayanagi prescription to compute en¬ 
tanglement entropy. The formula was hrst conjectured in [50, 51] and later established 
rigorously in [52], 

n The minus sign in (2.72) accounts for the fact that an increase in pressure generally induces a 
reduction in volume. 
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Let us start with some field theory generalities. At zero temperature, the system 
is in a pure state and can be described by a single state vector T), and by a density 
matrix p = |d/)(\l/|. The von Neumann entropy of the system is defined to be 


S = -tr[(5logp] . 


(3.1) 


We can consider a Cauchy surface in spacetime, E, and let A and A c be a subregion of 
E and its complement, respectively. The Hilbert space factorizes into H = Ha ® Ha c , 
and one may define a reduced density matrix for region A by tracing out the degrees of 
freedom of the complement, 

Pa = tr ac [p]. (3.2) 

The entanglement entropy of region A is then defined as the von Neumann entropy of 
the reduced density matrix, 

S A = -tr [p A log p A ] ■ (3.3) 

Consider now the finite temperature case. For a local quantum field theory, the system’s 
Hamiltonian H can be split into H = Ha + Ha = + Hqa, where Ha is the Hamiltonian of 
the degrees of freedom in the interior of A, Ha c is the Hamiltonian of the complement 
of A and Hqa is the Hamiltonian of the interface of A with its complement. The total 
density matrix is then given by: 

P = e~ pH , (3.4) 

and the definitions of the von Neumann entropy and entanglement entropy are the same 
as in the zero temperature case. 

Let us now move to the holographic prescription. According to the Ryu-Takayanagi 
formula, Sa is computed a la Bekenstein- Hawking, 


Sp, = 


Area(TA) 


4 G 


(3.5) 


N 


where Ta is a codimension-2 minimal surface with boundary condition <9T a = dA, and 
Gn is the gravitational Newton’s constant. 12 Notice the close similarity with equation 
(1.2). Indeed, if we consider A to approach the full E, it is possible to show that part of 
Ta tends to wrap the horizon of the black hole and the leading term of the entanglement 
entropy is dominated by thermal entropy, as expected. 

Finally, note that in a compact space, the complement of A basically has the same 
shape as A so one might naively think that the same minimal surface computes both Sa 
and Sa c - This is indeed true in the vacuum because in any pure state Sa = Sa =■ However, 
in a mixed state this is no longer true. In this case there is an additional topological 


12 The notion of a minimal surface only makes sense in Euclidean space. Therefore, the formula stated 
above only works for static spacetimes, where, by symmetry, the minimal surface lies on a constant time 
slice. For time-dependent backgrounds, the prescription can be made covariant by replacing the minimal 
surface condition by an extremal surface condition [53]. 
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requirement for the minimal surface in the bulk [54]: has to be homologous to A. If 

there is a horizon in the bulk, the minimal surface is homologous to A but not to A c . 
To satisfy the homology requirement for Sa c , we have to include a second component: 
the surface of the horizon itself. For more details on the homology requirement, see for 
example [55, 56]. 


3.2 Entanglement entropy across the phase transitions 

In this subsection, we track entanglement entropy across the various phase transitions 
of the STU black hole and we will argue that entanglement entropy detects the Van der 
Waals transition whenever it occurs on the gravity side. 13 

Consider a constant time slice, and take the region A to be a spherical cap on the 
boundary delimited by 9 < 9 0 . We will pick a small value for 9 0 because we want to filter 
out the thermal part of the entropy. Specifically, we will choose 9 0 = 0.15. Now, since 
entanglement entropy is UV-divergent (as is the area of the minimal surface), it needs to 
be regularized. There are several ways to accomplish this; we will do it by subtracting 
the area of the minimal surface in pure AdS. In order to do so, we first integrate the 
area functional to some cutoff 9 C < 9 0 . Then we turn off the mass and all the charges of 
the black hole to obtain pure AdS in global coordinates: 

ds 2 = — (1 + g 2 r 2 )dt 2 H-— + r 2 dfl 2 . (3.6) 

1 + g z r z 

The minimal surface which goes to 9 = 9 0 on the boundary is given by 


r(9) = 


1 

9 



(3.7) 


We can easily integrate the area functional of this minimal surface up to 9 C . Then we 
subtract this quantity from the black hole one, and send 9 C to 9 0 . This gives us the 
renormalized entanglement entropy, which we will denote by A Sa- In particular, for the 
numerical computations we will choose 9 C to be 0.149. 

The plots generated for the entanglement entropy are shown in Figure 5. Let us 
start the discussion with the AdS-RN case, which was treated in [36] but is included here 
for comparison’s sake. In Figure 5(d), we plot A Sa versus f3 for three different values of 
the charge (above, at, and below the critical charge) and confirm that there is indeed a 
Van der Waals transition. Moreover, the critical charge and critical temperature agree 
with those for the TS transition —see Figure 3— and with the formulas (2.40). Since we 
had already established that the TS transition captures the same physical information 
as the PV transition, we can conclude that the entanglement entropy also has access to 
the extended phase structure. For large enough regions this is indeed expected, since 

13 See [57-59] for related works on holographic entanglement entropy and phase transitions. 
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Figure 5. Plots of entanglement entropy versus /3. Panel (a): The single-charge case, with 
Q = 2 (green), Q = 1 (orange) and Q = 0.5 (blue). Panel (b): The two-equal-charge case, 
with Q = 2.5 > Q c (green), Q = 2.0704 = Q c (orange) and Q = 1 < Q c (blue). The critical 
temperature is /3 C = 44.4187. Panel (c): The three-equal-charge case, with Q = 1.5 > Q c 
(green), Q = 1.0478 = Q c (orange), and Q = 0.97 < Q c (blue). The critical temperature is 
/3 C = 38.8714. Panel (d): the AdS-RN case with Q = 0.9 > Q c (green), Q = Q c = 5/6 (orange) 
and to Q = 0.75 < Q c (blue). The critical temperature is j3 c = 107ry / 3/2. 


entanglement entropy is dominated by thermal entropy. However, the non-trivial result 
we point out here is that this behavior holds true even for arbitrary small regions. 
Indeed, we find that the critical charge and temperature are independent of the size 6q 
of the boundary disk. 

The results for the three-equal-charge case are shown in Figure 5(c). As in the AdS- 
RN case, the three-charge case also presents a Van der Waals transition, and the ASA 
versus (3 plot is qualitatively similar to Figure 3(c). The critical temperature and critical 
charge also agree with the ones obtained from Figure 3, providing further evidence on 
our previous conclusion. Unlike the two previous cases, the two-equal-charge case and 
the single-charge case do not exhibit the Van der Waals transition, again like in Figure 3. 
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Figure 6. Plot of log |/3 — (3 C \ versus log | A Sa — AS/yd based on the 32 points closest to the 
inflection point of Figure 5 (c), ranging from AS a = 0.065 (the critical value) to AS a = 0.08. 


These plots are shown in 5(b) and 5(a), respectively. For the two-equal-charge case there 
is an unstable branch at subcritical charge, in analogy with Figure 3(b); this branch is 
squeezed out when the critical charge and the critical temperature are attained, with 
values that agree with those in Figure 3. Finally, the single-charge case does not have a 
Van der Waals transition, like in Figure 3. 


3.3 Critical exponent of entanglement entropy 

Given that the entanglement entropy behaves very similarly to the black hole entropy, 
even for away from the thermal regime, it is natural to wonder if it is possible to define 
a “specific heat” at constant P for entanglement entropy: 


r T 

Cp - T \~w~ 


1 

P 


(3.8) 


as well as an associated critical exponent a. For the three-equal-charge case, we have 
plotted log |(3 — f3 c | versus log |A Sa — AS/yd- To produce this plot, we used the critical 
curve of Figure 5(c), and in particular the 32 points which are closest to the inflection 
point on this curve. A linear fit is shown in Figure 6 and yields: 


log |/1 - p c | = 2.91447 log | A S A - AS a , c \ + 9.665 (3.9) 

The slope is consistent with the value 3, and therefore the critical exponent a for entan¬ 
glement entropy agrees with that of the black hole entropy. In particular, it seems that 
this critical exponent is independent of the size 6$ of the entangling region. However, 
there could potentially be a weak dependence on 9 0 that can only be detected at higher 
numerical accuracy. We will leave this interesting issue for future work. 
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4 Conclusions 


The extension of black hole thermodynamics by uplifting the cosmological constant to a 
dynamical variable is a fascinating topic that has generated a lot of recent interest. Ac¬ 
cording to this framework, black holes seem to undergo critical phenomena very similar 
to ordinary matter. For instance, the extended phase structure of the AdS-RN black hole 
reinforces the well-known fact that this system exhibits a Van der Waals-like liquid-gas 
transition. Similarly, other AdS black holes are shown to exhibit a very rich structure 
that includes a whole catalogue of classical phase transitions, such as reentrant phase 
transitions and multicritical points, making them a useful tool for holographic applica¬ 
tions of condensed matter systems. 

In this paper, we have investigated the extended thermodynamics of the 4-dimensional 
STU black hole in the fixed charge ensemble. We uncovered an interesting PV phase 
structure that interpolates between the Hawking-Page transition of the Schwarzschild- 
AdS solution and the Van der Waals transition of AdS-RN, depending of the number of 
charges turned on. When only one charge is present, we still have a Hawking-Page-like 
transition. When two charges are turned on, we have an intermediate phase structure 
that resembles the Hawking-Page transition at sufficiently low temperatures. Finally, 
when three or four charges are turned on, we have a Van der Waals-like transition. One 
might wonder if the nature of these two transitions is the same, since they both have 
a qualitatively similar behavior. In order to answer this question we have extracted 
numerically two of the critical exponents that characterize the transition for the three- 
charge case, and confirmed that they are indeed the same as those of the AdS-RN black 
hole. Therefore, the three-charge and AdS-RN black holes seem to belong to the same 
universality class, although it would be interesting to extend our analysis to all the other 
critical exponents. We will leave this interesting question for future work. 

We also showed that the transitions appearing in the PV space are also present in 
the TS diagrams and explained it by arguing that the relevant parameter is always a 
dimensionless combination between the temperature and the pressure, Q = /3 2 P. Thus, 
we can approach the transition either by varying P or T and keeping the other quantity 
fixed. The interesting point here is that, in terms of the dual field theory, these two 
process have completely different physical meaning; in one case one is changing the field 
theory while in the other case one is changing the state. This is due to the fact that 
by varying the cosmological constant we are changing the rank of the gauge group of 
the boundary theory, N, and hence the PV diagram can be interpreted as a space of 
theories. We argued that the thermodynamical volume V has a natural interpretation as 
a cutoff scale A*, so that the isotherms in a PV diagram are interpreted as the running 
of N( A*) with the energy. In this context, then, the phase transitions are interpreted as 
a requirement for the monotonicity of N( A*) along the RG-flow. 

Now, given that thermal entropy is subject to the same phase transitions of the ex- 
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tended PV phase space, it is natural to wonder if the same phenomenon can be observed 
for other field theory observables. One natural possibility to consider is entanglement 
entropy, which is known to reduce to thermal entropy for large enough regions. In¬ 
deed, previous studies have shown that entanglement entropy also exhibits a Van der 
Waals-like transition for the AdS-RN black hole, but no connection with the extended 
thermodynamics was pointed out. Remarkably, we have shown that even for entangling 
regions of small size, the Van der Waals transition is still present, and it happens at 
exactly the same critical temperature, pressure and charge. This shows that entangle¬ 
ment entropy have also access to the extended thermodynamics and can be used as an 
efficient tool to diagnose its phase structure. In addition, the same statement seems 
to hold for other charge configurations of the STU black hole, regardless the nature of 
the transition. This provides further evidence supporting our results and makes the 
conclusion more robust. Finally, we computed numerically the critical exponent for the 
entanglement entropy and showed that it coincides with that of the thermal entropy 
whenever the Van der Waals transition is present. 


A number of possibilities for the extension of our work can be considered. First, 
we can perform a similar study in backgrounds that exhibit peculiar properties. One 
of such examples is the AdS-Taub-NUT geometry, which is known to have negative 
thermodynamical volume [26-28] . It would be interesting to see how this case would fit 
in the interpretation presented in this paper. Similarly, there are other examples that 
belong to this category. Two of them are the AdS-Kerr-bolt black hole and the Lifshitz 
black hole [60, 61]. In the first case, the issue arises because the volume is not the same 
as the one arising from naive geometrical considerations while, in the second, there is an 
ambiguity in expressing the cosmological constant lcngthscale in terms of the horizon 
radius lcngthscale. We can also consider the extended thermodynamics in dynamical 
settings, e.g. in Vaidya backgrounds, at least in the adiabatic limit. In such a regime 
one can still define an effective thermodynamics in terms of an apparent horizon and 
one might encounter interesting phenomena, specially in charged backgrounds [62-65]. 
A second line of study would be to investigate whether these phase transitions can be 
observed or not in other field theory observables. Some natural possibilities to explore 
are Wilson loops, correlation functions, or other entanglement-related quantities such as 
mutual information or renormalized entanglement entropy (a good starting point would 
be [66]). Finally, it would also be interesting to explore how 1/N corrections affect 
the phase structure of the extended thermodynamics and the entanglement entropy. 
Since our results are strictly valid in the large -N limit, one might wonder if the phase 
transitions we observed here can receive corrections, or even disappear (see e.g. [67]). 
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A Holographic heat engines 

Treating the cosmological constant as a dynamical variable grants us access to the PV 
phase space which, as explained in Section 2.4, can be understood as a space of field 
theories with different N and cutoff scale A*. Thus, we can think of constructing several 
cycles by moving in this RG-flow along different paths. Such processes can be thought 
of as holographic “heat engines”, in the sense of [35]. In this section, we present a few 
analytic results about STU black holes as heat engines, especially the efficiency of these 
systems in the high pressure regime. The results presented here are an extension of the 
paper [35], where the efficiency of the AdS-RN case was studied. 14 

For any cycle, the efficiency 7] is defined to be: 


W 


(A.l) 



where W is the work generated during one cycle, and Qh is the net heat flowing into 
the engine. Since W = Qh — Qc , where Qc is the waste heat leaving the engine, the 
efficiency can be equivalently defined as: 


In a reversible process, the entropy leaving the engine must equal the entropy flowing 
into the engine. Since Q = TS, this implies 


Qc — TfrQn , 

i-H 


(A.3) 


and therefore the efficiency of a reversible process takes the form 



(A.4) 


14 Other holographic heat engines were constructed recently in [68]. 
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This efficiency is called the Carnot efficiency, and is the maximal value that any heat 
engine can achieve. Here, we will verify that none of the STU black holes can violate 
the Carnot bound for the efficiency. 

Let us focus on the large-P limit or, equivalently, large -N limit. For the single- 
charge case, the high pressure expansions of the temperature, entropy and volume are 
(working to second order in g~ 2 )\ 


, r 3 2 ,2 Q 2 + r 2 + 3Q A + 2Q 2 r% , 

T — - — r + g H--—3-——-h • • ■ 

47T 47rr 4 2iTg 2 r\. 

2 , 27r( 5 2 1 ( 2ttQ 2 , lOvrQ 4 ^ 

S = nr + + ^ + —— ) + • • 


4 3 47 tQ- 1 f AttQ 2 167tQ 4 \ 
3 + r + g 2 g A \ r\ r + J 


+ 


r+g 2 g A V r + 

For the two-equal-charge black hole, the expansions are 

3 2 , 4 Q 2 + r\ 2(4Q 4 + Q 2 r 2 + ) , 

T — —r + g H-- —3 -- V ■ ■ 

47T 47rr^_ ng z r+ 

c_„ 2 , 47 r( 5 2 1 ( 4ttQ 2 , 32t tQ 4 \ 

S — nr + + + r a j + • • • 


r 2 n 2 

4 3 87 tQ 2 1 f 8nQ 2 1607t(5 4 \ 

= 3’ rr+ + 7^“?iV~ + _ 3d _ J 


r+g 2 g 4 V 


+ 


Finally, for the three-equal-charge black hole, the expansions are: 

3 2 6 Q 2 + r 2 + 3(13Q 4 + 2Q 2 r 2 + ) 

T ~ JZ r +9 + —7 —3 7^-7 + ' 


(A.5) 

(A. 6 ) 

(A.7) 

(A. 8 ) 

(A.9) 

(A. 10 ) 

(A.ll) 
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V 


+ 


+ 


(A- 12 ) 


(A.13) 


If we keep only the leading term in each of the expansions above, then we obtain an 
approximate equation of state valid in the high pressure, high temperature regime: 




1/3 


T 


3 J V 1 / 3 ' 


(A. 14) 


Comparing with the equation of state of the AdS-RN solution, we find that the right- 
hand side of (A. 14) is the first term in (2.35). Therefore, all four cases of the STU 
black hole have the same equation of state in the high pressure, high temperature limit. 
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Figure 7. The isotherm-isochore engine. The curvy lines are isotherms, and the straight lines 
are isochores. In terms of the boundary theory, an isochore can be interpreted as a process 
where N is varied while keeping the cutoff A* fixed. 


However, due to the subleading terms, the entropy and volume in general no longer 
determine each other, unlike the AdS-RN case. For the single-charge case we have: 

1/3 

*-+ ^ / . vn \ 

V = 


s 3/2 + ( ^ 

1 32 


3a/ 7r 

for the two-equal-charge case: 

v=^ r s 3 / 2 +(— 
3v^F V 9 

and for the three-equal-charge case: 

4 


Q 4 
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+ 


1/3 


Q 4 
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+ • • • , 


V = —s 3/2 + 
3 V7T 


37T 2 


1/3 


Q 4 

p2yb/3 


+ 


(A. 15) 


(A. 16) 


(A. 17) 


Notice that the leading term on the right-hand side of (A. 15)-(A. 17) is exactly the 
Reissner-Nordstrom relation (2.31). The second terms have a common form but they 
differ by a multiplicative numerical constant. 

A.l Isotherm-isochore cycle 

In this subsection, we compute the efficiency of Stirling cycles, i.e. cycles consisting of 
two isotherms and two isochores —see Figure 7. In the case of the AdS-RN, the fact 
that S and V are proportional to each other implies that the adiabatic curves coincide 
with the isochoric ones, and Stirling cycles coincide with Carnot cycles. Therefore, the 
Stirling cycle for AdS-RN operates at the optimal efficiency rjc■ To see this, consider 
the work Qh, Qc done along the upper and lower isotherms: 

Qh = T h AS^ 2 = T h (^0 7 (R 2 2/3 - R, 273 ), 


4 


(A.18) 













Qc 



(A.19) 


where the subscripts 1, 2 , 3, 4 refer to the points labeled on Figure 7. Since Vi = V 4 
and V 2 = V 3 , we deduce that the Carnot efficiency is attained. 

Next, consider the single-charge black hole. Using the expansion (A. 15), we now 
fold that the heat flows along the isotherms pick up some corrections: 

Qh = Qf + T„*f- (pFg -. (A. 20 ) 

<30 = 13™ + T c 4 (y2p2 - \fpl) • ( A ' 21) 

where Q™ and are given by (A.18) and (A.19), respectively. Next, we use the high 
pressure equation of state (A. 14) to rewrite the correction terms above: 


Qh =QT + $- 

J-H 

Qc = ($? + %■ 

J-C 


It follows that the efficiency is: 
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(A. 22 ) 
(A.23) 


(A.24) 


where we have used the fact that V\ = V 4 and V 2 = V 3 . Since Tc < Th , it follows 
from the above that the efficiency is less than the Carnot efficiency. The efficiency for 
the two-charge and three-charge cases can be obtained from the above by rescaling the 
common charge Q and its found to be less than the Carnot efficiency as well. 


A.2 Isobar-isochore cycle 

Next, consider a cycle composed of a rectangle on the PV diagram, i.e. two isochores 
and two isobars, as depicted in Figure 8. For this cycle, it is enough to truncate the 
high temperature expansions to first order in g~ 2 . To this order, all the three special 
cases (one-, two- and three-equal-charges) can be treated in the same time by defining 
an effective charge Q e ff as follows: 

Q\ff = Q\ + Ql + Ql + Ql ■ (A.25) 

In particular, Q e ff = Q f° r the single-charge case, Q e ff = \f^Q for the two-equal-charge 
case, and Q e ff = V3Q for the three-equal-charge case. To compute the heat flow along 
the isobars, we need the specific heat at constant pressure: 

° p= t (I?) p =2wrp+ 3 y Q ‘" + t - < a ' 26 > 
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Figure 8. The isobar-isochore engine. The vertical lines are isochores, and the horizontal 
lines are isobars. In terms of the boundary theory, an isochore can be interpreted as a process 
where N is varied while keeping the cutoff A* fixed. In an isobar process, on the other hand, 
A* is varied while N is kept fixed. 


The heat flow at the upper isobar is given by: 

rT 2 r r +,2 ut ^ 

Qh= / C P (Pi, T)dT = / Cp( gi ,r + )—dr+, (A.27) 

Jt, J r +, 1 ar + 

and integrates to: 

Qh = 2 - r 3 +1 ) + hr +,2 - r+,i) + 2 Q 2 eff (— -—. (A.28) 

Z Z ' \ r +,2 r +,l/ 

The work done during one cycle is simply the area of the rectangle inside the cycle: 


w = (V 2 - V 1 )(P 1 - P 4 ) 


Explicitly, 


w = i(9? - gl)(rl , 2 - <,) + §<£„ (1 - I) ( 


r+,2 r +}1 


The efficiency is then found to be: 

W 


V = 


Qh 


= |l-4Wl- r +,2 — r +,l 


9i 




(A.29) 


(A.30) 


(A.31) 
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As expected, the efficiency is again smaller than the Carnot efficiency. Moreover, the 
leading correction to the efficiency is charge-independent. If we convert g and r + to P 
and S, equation (A.31) becomes 



(A.32) 


This agrees with equation (20) in [35] for the AdS-RN case. 
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